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Abstract 
T— I ■ 

^ ' We address the issue of lepton flavour violation and neutrino masses in 

CN . the "fat-brane" paradigm, where flavour changing processes are suppressed by 

! localising different fermion field wave-functions at different positions (in the 

' extra dimensions) in a thick brane. We study the consequences of suppressing 

■ lepton number violating charged lepton decays within this scenario for lepton 

, masses and mixing angles. In particular, we find that charged lepton mass 

I matrices are constrained to be quasi-diagonal. We further consider whether 

Q_i' the same paradigm can be used to naturally explain small Dirac neutrino 

^ . masses by considering the existence of three right-handed neutrinos in the 

Q I brane, and discuss the requirements to obtain phenomenologically viable neu- 

' trino masses and mixing angles. Finally, we examine models where neutrinos 

^> . obtain a small Majorana mass by breaking lepton number in a far away brane 

^ , and show that, if the fat-brane paradigm is the solution to the absence of 

^ ' lepton number violating charged lepton decays, such models predict, in the 

I absence of flavour symmetries, that charged lepton flavour violation will be 
observed in the next round of rare muon/tau decay experiments. 



* gabriela. barcnboim@cern . ch 

tgbranco@thwgs.cerii.ch and gbranco@cfif.ist.utl.pt 
^ degouvea@mail.cern.ch 

Wrebelo@thwgs.cern.ch and rebelo@cfif.ist.utl.pt 

^On leave of absence from Departamento de Ffsica, Instituto Superior Tecnico 
Av. Rovisco Pais, P-1049-001, Lisboa, Portugal. 



1 Introduction 



A few years ago, a new paradigm for addressing the gauge hierarchy problem was 
proposed In the so-called "large extra dimensions scenario," (ADD scenario) the 
explanation for the large discrepancy between the weak scale and the Planck scale 
is that the fundamental scale of gravity is actually close to the weak scale, and that 
gravity seems weak to us because it propagates in more, compact, dimensions, while 
Standard Model fields are bound to four-dimensional subspaces embedded into the 
higher dimensional world. 

One of the challenges for the ADD scenario is solving the well known problems 
of proton stability, absence of flavour changing neutral current processes, etc. The 
situation is particularly problematic because one has to consider nonrenormalizable 
operators which are only suppressed by powers of the "quantum gravity scale" and 
whose coefficients are completely unknown since they parametrise our ignorance 
concerning the physics above the ultraviolet cutoff. It is important to ask whether 
there are "low energy" mechanisms guaranteeing that most of these dangerous terms 
are appropriately suppressed, independently of the unknown beyond-quantum-field- 
theory physics. 

A novel solution to the higher-dimensional fiavour problems was proposed by 
Arkani-Hamed and Schmaltz (AS scenario). They argue that if we live in a 
"fat" four-dimensional spacetime (fat brane), it is possible to effectively suppress 
unwanted higher dimensional operators by localising different fermion fields to dif- 
ferent coordinates in the extra dimensions. As was also pointed out in 0, and 
further explored in [0, this solution also provides an interesting explanation for 
the hierarchy of the quark masses and the values of the CKM mixing angles and 
CP-violating phase. It should be pointed that the AS scenario leads to startling 
consequences P|, which may be observed at next-generation collider experiments. 

In this paper, we address the issue of lepton fiavour violating phenomena in 
charged lepton decays and neutrino masses in the AS scenario. The AS scenario is 
natural for solving fiavour changing problems, and indeed we find that it can very 
efficiently suppress charged lepton fiavour violating decays such as /x, r — I'j, III. 
One of the issues we explore here is whether suppressing such processes imposes any 
interesting constraint on lepton masses and mixing angles. We find, for example, 
that the charged lepton mass matrix is forced to be quasi- diagonal, and that almost 
all the mixing in the lepton sector must come from the neutral sector. Such a 
situation is certainly not observed in the quark sector (for an example, see [Q), and 
seems to already constrain some extra-dimensional neutrino mass scenarios. 

The AS scenario also provides a natural mechanism for obtaining small Dirac 
neutrino masses, which has not yet been explored in the literature. It is important, 
therefore, to study whether similar success is obtained in explaining, for example, the 
large leptonic mixing angle which is required by neutrino oscillation solutions to the 
atmospheric neutrino puzzle 0. We will show that in order to obtain large mixing 
angles, separations between different leptonic fields have to be closely related. The 
situation appears to be rather finely tuned when the fiavour constraints mentioned 
above are also imposed. 
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The paper is organised as follows. In Sec. 2, we address the problem of flavour 
changing neutral currents in the ADD paradigm and how it is solved in the AS 
scenario. We concentrate on flavour violating charged lepton decays. In Sec. 3, we 
address the issue of lepton masses and mixing angles in a straightforward extension 
of the AS scenario, namely, adding right-handed neutrinos to the fat brane. We first 
discuss how "naturally" the AS scenario can explain the large hierarchies observed 
in the quark and fermion sectors, and what is required in order to obtain large 
mixing angles. We then analyse two-family models in great detail, concentrating on 
the constraints that must be met in order to obtain large mixing angles. Finally, 
we discuss solutions to the neutrino puzzles in the case of three generations. In 
Sec. 4, we discuss a mechanism for obtaining small neutrino masses without right- 
handed neutrinos. Here, constraints from the charged lepton sector severely restrict 
solutions to the neutrino puzzles. Our conclusions are presented in Sec. 5. 

2 Lepton Flavour Violation and Large Extra Di- 
mensions 

It is a general principle of effective field theories that all operators consistent with the 
exact symmetries of Nature are present, including nonrenormalizable, "irrelevant" 
operators. These operators are understood as being suppressed by the large energy 
scale above which the effective field theory is no longer applicable. These include, 
for example, {1/ K^)QQQL^ which mediates proton decay and {1/A)LHLH, which 
breaks lepton number and gives the neutrinos a Majorana mass after electroweak 
symmetry breaking. It is also well known that A is constrained to be at least of 
the order of the grand unification scale {Mqut — 10^^ GeV) by the fact that the 
proton lifetime is larger than ~ 3 x 10'^^ years [0]. On the other hand, values of 
A ~ 10^~^^ GeV are required in order to correctly account for the atmospheric 
neutrino puzzle via neutrino oscillations. 

In theories where the "quantum gravity" scale is small (a few TeV), one is faced 
with the challenge of explaining why, even though A ~ 1 TeV, phenomena which 
violate the conservation of "accidental" global symmetries have not been copiously 
observed. The two irrelevant operators mentioned above can be eliminated by 
assuming, for example, that (gauged) B — Lis a. true symmetry of Nature (probably 
weakly broken, if one wants to explain the matter-antimatter asymmetry of the 
universe), but other dangerous operators still remain. These include operators 
that violate individual flavour number conservation and mediate phenomena such 
as meson-antimeson mixing and flavour changing neutral currents. 

One intriguing solution to these problems was proposed by Arkani-Hamed and 
Schmaltz 0. They postulate that the standard model fields are constrained to a 
"fat brane," which is infinite in three space directions and occupies a finite volume 
in n extra, orthogonal, compact dimensions. It is also postulated that gauge fields 

*We use the standard notation where Q and L refer, respectively, to quark and lepton SU{2) 
doublets, while U, D, and E denote, respectively, up-type antiquark, down-type antiquark, and 
antilepton SU{2) singlets. H is the SU{2) Higgs boson doublet. 
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and the Higgs scalar can freely propagate in the entire volume of the fat brane, 
but fermions are confined to specific "points." Finally, if one assumes that different 
fermion fields are located at different points in the extra dimensions, it is easy 
to see that some operators are "forbidden" by locality, e.g., an operator such as 
{1/ K^)QQQL can be very efficiently suppressed simply because Q's and Us live in 
different worlds, and don't "see" one another 0. 

In this section we discuss, within the scenario described above (AS scenario), 
what requirements are imposed on the distances between different leptonic fields in 
order to account for the (negative) experimental searches for charged lepton decays 
which violate lepton flavour number. Here we will concentrate on the lepton flavour 
violating muon decays 

IX — >■ 67, jjL —>■ eee, (2.1) 

and tau decaysjj] 

T /i7, T — > 67, r — > /i/x/i, r — > eee. (2.2) 
The process l'^ —>■ /'7 {k, I = t, fi, e) is mediated by the effective Lagrangian 

- U^i^, = ^ {i1^,AF^') + ^ {f^a.J.F^'') + h.c, (2.3) 

after electroweak symmetry breaking. Here the subscript L, R refers to the chirality, 
F'^'^ is the electromagnetic field strength, v is the Higgs vacuum expectation value, 
and a/j L are dimensionless couplings. 

The process I'' iHH^ is mediated by the effective Lagrangian (after Fiertz 
rearrangements) 

+^ + &2(/14)(44)+ 

+h{lhAMl%) + hHilAWRl'lR) + h.c] , (2.4) 

where 6's are dimensionless couplings. Eq. ( p.3| ) mediates l'^ iHH^ by attaching an 
iH'- pair to the photon leg. 

In the AS scenario, one can relate the constants a, h to the distance between 
different leptonic fields and to higher dimensional couplings by integrating out the 
extra dimensions. One obtains, for example, 

aROcanJ dye-'^'^y-^'^^' e-^^'^y-"^' oc a Re' -2'''^^''-"^' (2.5) 

and 

64 oc /34 / dye-'^'^y-'^^" (e-^'^^^-'')')' oc (3,e--^^"^'''-'^^\ (2.6) 



^We do not include the decays r — > fiee, e^/z in our discussion. The reason for this is that the 
constraints imposed by these channels are far weaker than those imposed by — s- ej, eee, as will 
become clear shortly. 



3 



assuming one extra dimension and defining the higher- dimensional wave-function 
of the fermions to be Gaussians of width /i~^/-\/2 centred at position Ik^i and ek,i 
for l£ and , respectively. The couplings a, (3 are dimensionless numbers, which 
are assumed to be of order one, since there is, a priori no reason for them to be 
either very small or very large. All a and (3 will be, henceforth, set to unity unless 
otherwise noted. 

From Eqs. ( p.3|J2.4| ) and assuming the AS scenario, it is easy to compute the 
branching ratios for the various lepton decays (see, for example, as a function 
of A (the "quantum gravity" scale) and the distance between different lepton fields. 

Fig. m depicts the branching ratio for /i — 67 (right) and r l^y, I = e, fj, 
(left) normalised by the branching ratio for /i evv and r — > Ivv respectively, as 
a function of AL. Two different values of the cutoff, A = 10 TeV and 1 TeV, are 
shown (top and bottom, respectively). Here, AL is either \l^ — ee\ or — le\ for 
the muon decay and |/t — q| or \er — li\ for the tau decay, in units of yU~^. 

Fig. ^ also depicts the branching ratio for yU — > eee (right) and r — > I = e, fi 
(left) normalised by the branching ratio for fi — * euu and r —>■ Ivv respectively, as 
a function of AL. Here, only nonzero 63 and 64 coefficients were considered. The 
reason for this is simple: given the requirements imposed by the non-observation of 
/i — >• 67 and r Z7, the contribution of a^^L, ^1,2,5,6 are strongly suppressed, while 
63,4 are, a priori, unconstrained. Therefore, in this case, AL is either \l^ — le\ or 
|e^ — Cel for the muon decay and 1/,- — ^zl or le,- — e^l for the tau decay, in units of 

Some comments are in order. First of all, it must be said that the AS scenario, as 
expected, is very efficient when it comes to suppressing rare lepton decays. Even for 
A = 1 TeV, the most stringent requirement (from — > 67) is that |/^ — Ce], |/e — ~ 
6.5/i~^ (proton decay constraints require distances between quarks and lepton of 
order 10/i~^ 0). Second, one may wonder how strict the constraints are since, after 
all, there are order one coefficients floating about. It is easy to see that things do not 
change significantly if these coefficients are allowed to vary. For example, varying 

from 0.3 to 3, the lower bound on |/r — from r — /i7 and A = 10 TeV changes 
by 17%. A remaining question is what do these bounds imply for lepton masses and 
mixing angles. This will be addressed in the next sections. 

Before proceeding, however, it is worthwhile mentioning at this point that other 
lepton flavour violating decays, such as K — >• /le, hadronic tau decays and /i-e 
conversion in nuclei^ have not been considered, since they can be safely suppressed 
by assuming that quarks and leptons are well separated (which is generically required 
in order to suppress proton decay), and that the distance constraints obtained in 
these cases do not play a significant role when it comes to addressing lepton masses 
and mixing angles. 



■I-The Lagrangian Eq. (2.4) mediates /i-e conversion in nuclei at higher order in QED (in the 
case of the magnetic moment operators) and at higher loop level and higher order in QED in the 
case of the four-fermion operators. The constraints imposed on the operators are not as strict as 
the ones from lepton flavour violating charged lepton decays. 
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Figure 1: Branching ratio for ^ ^ X (right) and r ^ X, (left) normahsed by the 
branching ratio for /i evv and r Ivv respectively, as a function of AL, for two 
values of the quantum-gravity scale, A = 10 GeV (top) and 1 GeV (bottom). See 
text for details. The horizontal lines indicate the current experimental constraints 
from searches for rare muon and tau decays 0]. In the case of tau decays, only the 
most and least stringent bounds are shown. 
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3 Right-Handed Neutrinos in the Brane 



Another intriguing possibility is using the AS scenario to explain the current hi- 
erarchy of fermion masses and mixing angles 0. Indeed, it has been shown that 
localising different quark fields to different points in the extra dimensions is an 
efficient and "natural" way of explaining not only all quark masses and mixing 
angles 0, but also the observed CP violation in the quark sector [Q. 

In this section, we study whether the same idea can be used to explain lepton 
masses. Charged lepton masses were already discussed in the context of massless 
neutrinos [Q, but no similar study has been performed for neutrino masses. 

The idea is quite simple. We assume the existence of three standard model 
singlet fermions, which we refer to as right-handed neutrinos (N), and assume that 
these are also localised at different points in the extra dimensions. If this is the 
case, neutrinos acquire ordinary Dirac masses after electroweak symmetry breaking, 
exactly like the quarks and charged leptons. In the AS scenario, the explanation for 
the smallness of the neutrino masses is straightforward: neutrino masses are tiny 
because A^'s and L's are separated by distances which are only slightly larger (as 
will become clear later) than, for example, the Q's and D's.|^ It will turn out, as 
we will discuss in detail later, that the most serious challenge for the AS scenario 
with right-handed neutrinos in the brane will be to naturally account for the large 
mixing angles in the lepton sector which are required in order to solve the neutrino 
puzzles p. 

Before proceeding, it is worthwhile to digress a little on the philosophy behind 
the AS scenario when it comes to explaining fermion masses and mixing angles. 
Perhaps the most attractive, and certainly the novel idea behind the AS scenario 
is the possibility of solving the flavour problem without postulating new, horizontal 
symmetries. Instead, it is simply assumed that the fermions are separated by 
distances of a few (in units of the width of the fermionic higher dimensional wave- 
function which could be chosen at random. All the work is done by the 
Gaussian suppression of the effective four dimensional Yukawa coupling (as discussed 
in the previous section. See also [Q, ^): 

^o,e-iM^(/.-/.)^^ (3.1) 

where fi are the positions of the different fermion fields {Q, U, D, L, E, N). 

In order to further quantify the "naturalness" concept, let us assume that, say, 
a Q and a U fields are localised in one extra dimension, of size L = 20/i~^ (see 
for a discussion on the expected size of Lfi). If the positions are chosen completely 
at random]] it is easy to show that the probability of finding \q — u\fi between some 

* Another possibility for obtaining naturally small Dirac neutrino masses is to assume that the 
neutrinos couple to right-handed neutrinos which live in the bulk Lately, this possibility has 
received lots of attention and will not be discussed in this paper. 

^It has recently been discussed in that the relative distance between different fermion fields 
may be significantly altered by gauge interactions. We will disregard this effect for the sake of 
simplicity. 
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Figure 2: The probability for obtaining a suppression of the higher dimensional 
Yukawa coupling of a given order of magnitude. See text for details. 

value 6 and 5 + d5 is given by P{6)d6, where 

P{S) = 0.1 - .005(5. (3.2) 

It is straightforward to compute the probability that the suppression due to the 
small overlap between the two extra dimensional wave-functions is of a given order of 
magnitude. Fig. ^ depicts the probability that the suppression factor is in some order 
of magnitude bin. It should be noted that while the figure only depicts suppression 
factors larger than 10~^°, the probability that the suppression factor is smaller than 
10"^° is almost 10%! 

As expected, very tiny suppression factors can be obtained, even if all fields are 
constrained to lie relatively close to one another. An important fact is that, while 
random choices for the positions (in one extra dimension) seem to favour effective 
Yukawa coupling which are suppressed by an order of magnitude with respect to 
their higher dimensional counterparts, the probability of obtaining a large spread of 
suppression factors is quite large. It is this property which renders the AS solution 
to the fermion mass hierarchy attractive - random choices for the distances seem to 
provide a large spread for values of the Yukawa couplings, and therefore the observed 
hierarchies are naturally explained. 

One can further pursue this logic and see what are the consequences for the 
mixing angles. We use a toy model of only two families of up-type quarks as an 
example, but it will become clear in due time that it is easily generalised to include 
three families and down-type quarks. 
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The mass matrix can be written as 

M, = p[l (3.3) 

where p is roughly the top quark mass and a,b,c,d < 1 are the suppression fac- 
tors. The message from Fig. ^ is that, "normally" one expects a ^ 6 3> c S> o? 
(for example). If this is the case, one can trivially show that the matrix that 
diagonalizes MuM}j is parametrised by one mixing angle, 6', such that cos 26' = 
±1 =F 0{(? / a?' ,cbd/ a^), i.e., 9 is either close to or 7r/2, such that the mixing is 
quite suppressed. 

While this situation fits very nicely the quark masses and mixing angles, it 
poses serious problems for leptons, given that, according to the neutrino oscillation 
solution to the neutrino puzzles, at least one of the leptonic mixing angles is "large." 
This is in sharp contrast to standard four dimensional approaches to the fermion 
mass puzzles. There, quark masses require a significant amount of underlying 
structure in order for one to understand the hierarchy of masses and the small 
mixing angles, while it has been argued that neutrino masses can be obtained from 



operators with random order one couplings |jT3| . The same is true in the case of the 



"simplest" neutrino mass model [T^, where neutrino masses are obtained via the 
introduction of one single right-handed neutrino and the use of l/Mpianck operators. 

In order to obtain large mixing angles, it is imperative that some elements of 
the mass matrices are of the same order. Using Eq. ( p.3|) as an example, it is again 
simple to show that, if one wants 0^7r/4, a~c (for example) is required. 

The challenge is, therefore, to obtain a mechanism satisfying this condition. 
Flavour symmetries, for example, would certainly do the job, but they would spoil 
the most attractive feature of the AS scenario. A different option, which still 
maintains the extra-dimensional nature of the AS scenario, is to assume that, as 
a consequence of the underlying localising mechanism, some fields are positioned on 
top of each other, and that some distances are, therefore, "naturally" very similar. 

With this in mind, we proceed to discuss how realistic lepton masses and mixing 
angles can be obtained in the AS scenario. We first discuss the simpler two family 
case, and further address the issue of how "finely tuned" certain parameters have to 
be in order to obtain large mixing. Later, we discuss solutions for the more involved 
case of three families. 



3.1 Two Neutrino Families 

In this subsection we address in detail the case of only the two "heavy" lepton fam- 
ilies, concentrating on the observed large mixing which seems to be required in the 
p — T sector according to neutrino oscillation solutions to the atmospheric neutrino 
puzzle. It proves convenient to consider separately the different configurations which 
potentially yield maximal mixing in order to understand what constraints over the 
space of locations have to be met in order to obtain the appropriate mass-squared 
difference and mixing angle. These configurations are depicted in Fig. |^. 

Two important comments concerning our discussion are in order. First, we 
assume that the order one coefficients which accompany the exponential suppressions 
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Figure 3: Different two family configurations which potentially lead to maximal 
mixing. 

(the higher dimensional Yukawa couplings) play no significant role, and set all of 
them to unity. We comment on the relevancy of this approximation later. Second, we 
impose that there is no significant mixing coming from the charged lepton sector. It 
turns out that, specially when considering only the fi — r sector, the non-observation 
of r — i> /i7 combined with the requirement that the obtained r mass agrees with the 
experimental value automatically forces the mixing in the charged lepton sector to 
be small. Note that the constraint of small mixing in the charged lepton sector will 
be considered a posteriori, and its consequences analysed for each of the different 
configurations. 

By imposing the t ^ fi'-f constraint as discussed in the previous section, one 
obtains - e^|, |/^ - e^| > 3.73/1-^ for A = 10 TeV (5 TeV), which 

in turn means that the off-diagonal elements of the (two generations) Dirac mass 
matrix for the charged leptons are less than (approximately) pe~^'^^^^ = 88 MeV 
for p = 1.5 X rrit, where = 166 GeV is the top quark mass (for details concerning 
the choice of p, see [Q). This, in turn, implies that one of the diagonal elements is 
necessarily of the order of the tau mass, while the other elements are at most of the 
order of the muon mass. In summary, 

fe(l + bce) be \ ,^ 

m^-^m^[^ C6 l-f(62 + c^)j' (3.4) 

where b, c are numbers which are at most order one, e = m^/mr-, and m^^^ are the 
muon and tau masses, respectively. Note that Eq. ( 3.4|) is written in such a way 
that its eigenvalues are and m,- up to 0{e^). In this case, it is easy to estimate 
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the mixing angle, namely, 



sm 



e ~ < 10- 



(3.5) 



which is quite small. 

A fact that will prove very important in the future is that, even though there 
is no direct significant limit on — /t| from r — ^ /i/i/i decays (see Fig. [^), it is 
geometrically constrained because \It — er\ < — er\. It is easy to show that 
\ifj. — irl > — Gr \ — \It — Gr\ — (4.0 — 3.14)/i~^J^ independent of the number of 
extra dimensions. Therefore, as a guarantee that there are no large flavour changing 
effects in the charged lepton sector, a good rule of thumb is to keep — /t| ~ 
in which case the mixing in the charged lepton sector is necessarily small. 

Case A - Fig. §(A) leads to the neutrino Dirac mass matrix 



(3.6) 



where / = I2 and S = h — h- We restrict ourselves to 5 > 0, without loss of generality. 
As mentioned before, we will fix p = 1.5 x 166 GeV. 

It should be noted that the interchange of with rir (see Fig. |^) corresponds to 
a weak basis transformation which interchanges the two columns of niy but leaves 
rriuml unchanged. Note also that the configuration depicted in Fig. |5(A) leads 
automatically to very small mixing in the charged lepton sector, once the appropriate 
sizes for the neutrino masses are imposed. 

Diagonalizing m^ml starting from Eq. (p.6|) one obtains the mixing angle 9, 
which, to a very good approximation (assuming 5, e <^ /), is given by 



tan 26* 



where 



D2 



exp 

exp 
exp 
exp 



2N^{1 + N2) 
iD,-l)iD2 + l] 



2 X 

y(2/(5 + (5') 
y(2/e + e2) 



(3.7) 



(3.8) 

(3.9) 
(3.10) 
(3.11) 



The mass squared difference is 



Am^ = 2p V^''' ATi (1 + A^a) W 1 + 



tan2 26 ' 



(3.12) 



■fNote that mr — 1777 MeV corresponds to a distance of 3.14/i and = 105.7 MeV to a 

distance of 
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Table 1: Case A - The allowed range of /le and the maximum allowed value for 5/1, 
for different values of yu/, which satisfy 1 x 10"^ eV^ < Am^ < 8 x lO^^ eV^ and 



tan 2^1 > 2.3 15 



/i/ 


/ie 




7.1 


~ 1.1 - 1.3 


10-6 


7.2 


~ 0.9 - 1.2 


10-^ 


7.3 


~ 0.7- 1.0 


10-4 


7.4 


~ 0.5 - 0.8 


10-4 


7.5 


~ 0.3 - 0.6 


10-3 


7.6 


~ 0.1 - 0.4 


10-3 


7.7 


~ 0.0 - 0.2 


10-3 



Imposing large mixing, one can safely write 

Am^ ~ 2p2e-^'''A^i (1 + N2) . 



(3.13) 



Note that in this approximation Am^ is equal to the numerator of tan 2^ times 
p^e"^ ' . This important observation is valid in all cases we will discuss henceforth, 
as long as tan^ 26' ^ 1 (which is the condition we would like to meet). 

We start by obtaining bounds for /. Since the largest element of my is the (22) 
element, imposing an upper bound on the heaviest neutrino mass of ^ 10 eV^ 
leads to / > 7.08/i~^ On the other hand, from Eq. ( p.l3| ) and noting that N] , N', 



2<1, 



one obtains a rough upper bound for I by requiring Am? > 10-3 eV^: / ^ 7.8/i"^. 
This can also be derived directly from Eq. ( p.6| ) by noting that the largest neutrino 
mass matrix element should be (roughly) larger than the lower bound on y/Am? 
imposed by the atmospheric neutrino data. It is important to comment that for 
values of fil close to the upper bound mentioned above, the neutrino masses are 
hierarchical, while for fil much smaller than 7.7, the eigenvalues of the neutrino 
mass matrix are degenerate and much larger than V Aw?. 

Next, we determine what values of e and 5, for fixed /, are required in order 
to obtain large mixing and the appropriate Am? range. For fixed /, Am? fixes the 
order of magnitude of 2A'^i(l + A'^2) (see Eq. (|3.13| )), which is the numerator of tan 2^ 
(see Eq. ( p.7D ). Large mixing requirements then determine the order of magnitude of 
{Di — 1) (note that D2 + 1 is 0(1)), which allows one to compute an upper bound for 
5. Once this is done, one revisits the previous expressions and computes e. Table |I| 
contains values of /xe and {6/ 1) max as a function of fil. 

In order to examine how "finely-tuned" this scenario is, we look at the maximum 
value for the ratio 6/1, which is how much the distance \l^ — nfj^\ is allowed to deviate 
from \l^ — nT-\, normalised by the smallest distance. One can see that (5//)max ~ 10-3 
over all possible values of l,e, for small e. It is curious that the fine-tuning gets 



§For concreteness we consider large mixing to be | tan20| > 2.3, which agrees with the 99% CL 
contours obtained by the neutrino oscillation analysis of the SuperKamiokande atmospheric data 
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more severe as e increases. This can be understood easily by the following line of 
reasoning: /S.rn? essentially constrains e due to the fact that, it turns out, e ^ 5; 
for small values of /, a large suppression is required in order to obtain the correct 
value for Am^, which must come from roughly Hence bigger values of e are 

required. However, larger values of e suppress the numerator of Eq. (p.7|), requiring 
a larger suppression in the denominator, which is mostly a function of / and 5, and 
the largest allowed value for 5 decreases. 

Case B - Fig. ^(B) leads to the Dirac mass matrix 

where I21 = I and I22 = I + 5. We first assume 5 > 0. This leads, to a very good 
approximation (for 5 <^ Z), to 



2exp -^{2V + 215 + 5^ 

tan 2^ 



,2 



exp[— yU^Z^] — exp[— /i^a^] — exp[— /i2(Z + 5)^] 
and, remembering that we will require tan^ 29 to be large. 



(3.15) 



Am^ ~ 2p^ exp 



■ ^{2f + 215 + 5^) 
2 



(3.16) 



It is clear that the same approximate bounds on /, which applied in Case A, 
also apply here, and again we analyse what constraints are imposed on the other 
parameters for fixed values of I within the allowed range. It should also be noted 
that, similarly to Case A, the charged lepton mass matrix is safely diagonal, given 
that, it turns out, Z215Z22 ~ 7/i~^. 

For fixed Z, the Am^ constraint fixes the range of 5 and the numerator of tan 2^ 
(Eq. ( |3.15| )). The smallness of this term will determine how much the terms in the 
denominator of tan 2^ are required to cancel in order to obtain large mixing. The 
only free parameter which remains is a. Table |^ illustrates the severe fine-tuning 
between the value of / and a which is required in order to obtain large mixing. 
Only for / close to its upper bound is 5 small enough to allow for a near perfect 
cancellation between the terms ' and in Eq. (|3.15|) , in which case a 

plays no significant role, provided that ^ is small enough (which is satisfied for 
lia ^ 7.8). In this case, therefore, there is no need to tune the values of a and I22 
(this happens for /i/ = 7.7 and ^5 ^ 0.05). This does not mean, however, that there 
is no fine-tuning, since one is required to tune the values of /21 and I22 (for yuZ = 7.7, 

(5/Omax ^ 10-=^). 

Next, we examine the case Z21 > ^22? ^-C- ^ < 0. In this case, the dominant 
term in the denominator of Eq. (|3.15| ) is ('-I'^D ^ which cannot be cancelled by 
appropriately choosing a. The only way to obtain large mixing is, therefore, to 
choose <^ / and a large. Numerically, we obtain that cannot be larger than 
0.05. The bounds on both Am^ and tan 2^ can only be simultaneously satisfied for 
/i/22 = 7.6, 0.04 < n\5\ ^ 0.05 and /i/22 = 7.7 and ^ ji\5\ ^ 0.05, while /la ^ 7.8. 
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Table 2: Case B - The allowed range of jid and the maximum allowed value for 
\a — l\/l, for different values of /i/, which satisfy 1 x 10~^ eV^ < Am^ < 8 x 10"^ eV^ 
and I tan 2^ I > 2.3 H 



/i/ 




( a - / //)max 


7.1 


~ 1.0 - 1.3 


10-*^ 


7.2 


~0.8 - 1.1 


10^^ 


7.3 


~0.6 - 0.9 


10-5 


7.4 


~0.4 - 0.7 


10-^ 


7.5 


-0.2 - 0.5 


10-3 


7.6 


-0.05 - 0.3 


10-2 


7.7 


-0.5 - 0.1 


10-2 




-0 - 0.05 


unconstrained 



In this case, the fine-tuning occurs between Z21 and Z22, which have to be equal to 
one another to more than one part in 10^. 

Case C - Fig. |^(C) leads to the Dirac mass matrix 



rrty 



P 




which yields the approximate relations (e, 6 <^ I) 

2N[[1 + N^] 



tan 29 



and 
where 



{l-D[){l-D„ 



K 
D' 



exp 

exp 
exp 
exp 



-1?{215 + 5^ + 5e) 
-fi\2l6 + 6^ 
-l?{2le + e^] 



(3.17) 



(3.18) 
(3.19) 

(3.20) 

(3.21) 
(3.22) 
(3.23) 



Note that the terms in the square bracket in Eqs. ( p.l8| , |3n9D are, by definition, 
between 1 and 2 and play no significant role in order of magnitude discussions. As 
in the two previous cases, / is bound between roughly 7.1fi~^ and 7.8/^-^ (assuming 
the "atmospheric" Am^ and that ml < 10 eV^) and again we discuss the required 
cancellation among parameters for fixed values of I. As before, the requirement that 
Am^ satisfies the atmospheric neutrino data determines the order of magnitude 
of the numerator of tan 26 and, in this case, also determines the value of e. The 
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Table 3: Case C - The allowed range of yue and the maximum allowed value for 5/1, 



for different values of /i/, which satisfy 
I tan 2^1 > 2.3 ^ 



1 X 10-^ eV^ < Am^ < 8 x 10"^ eV^ and 
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fll 


fie 


(*^/Oinax 


7.1 


~ 1.1 - 1.3 


10-6 


7.2 


~ 0.9 - 1.2 


10"^ 


7.3 


~ 0.7- 1.0 


10-^ 


7.4 


~ 0.5 - 0.8 


10-^ 


7.5 


~ 0.3 - 0.6 


10-3 


7.6 


~ 0.1 - 0.4 


10-3 


7.7 


~ 0.05 - 0.2 


10-1 




~ - 0.05 


unconstrained 



requirement |tan26'| > 2.3 is met by choosing (1 — D[) very small, which implies 
that 6 has to be very small. Table ^ contains the value of /iZ, fie, and (5//)max which 
are obtained when trying to fit large mixing in the 2-3 sector with Am^ ~ 10"^ eV^. 

It is noteworthy that in the limit e — there is no constraint on 6. In this limit, 
the two lepton doublets are on top of each other, and large mixing is guaranteed, 
regardless of the location of the singlet fields. This situation, however, is in sharp 
contradiction with the bound from the flavour changing r — /i7 decay, as discussed 
earlier in this section. Indeed, when we fit for the charged lepton masses imposing 
that r — > /i7 is small, we find that fie ^ 1, and large mixing with the appropriate 
mass squared is only obtained for fil < 7.3 (and (5//)max ~ 10"^!). 

Case D - Fig. ^D) leads to the Dirac mass matrix 



2 

e 2 



(3.24) 



which leads to the relations 
tan 26* 



2N';[l + N'^] 



{D'i - l){D'i + 1) - (1 - e-^f'^^P)D'iD'i'- 



(3.25) 



and (for tan^ 29 large) 



where 



Am2^2pV^ ' Ar('[l + A/^'], 



(3.26) 



exp 

exp 
exp 
exp 



-fi\2ip + + aP) 
-fi^{2la + 
-fi\2l(3 + (3^\ 



(3.27) 

(3.28) 
(3.29) 
(3.30) 
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As in the previous cases, fil is naively bound between ~ [7.1, 7.8] simply by requiring 
the largest neutrino mass squared to be less than 10 eV^ and Am^ > 10~^ eV^. It 
is again useful to constrain the other parameters for fixed value of fil. In order to 
obtain a correct value for Am^, each value of / requires a specific range of a (which 
determines N", see Eq. ( p.26| )). Once a is fixed, we try to choose /? such that the 



mixing is large. In this case, however, it is easy to see that changing P for fixed a 
and / does very little to increase tan 2^, as it mainly affects the value of D2, which 
plays a diminished role when it comes to determining the order of magnitude of 
tan 2^ (see Eq. (g]^). 

The only possibility for obtaining large mixing is to choose a such that (-D" — 
1) is very small, i.e to choose a very small. Note that this choice automatically 
renders the term in the denominator of Eq. ( |3.25| ) proportional to 1 — e^^ "'^ small.|^ 
Small values of a require / to be large enough so that an appropriate value of 
Am^ is obtained. As a result, only values fil ^ 7.6, /la ^ 0.05 potentially yield 
phenomenologically allowed values for Am^ and tan 2^. 

As in Case C, the charged lepton sector imposes (here very severe) constraints 
over the allowed values for the — lr\ distance, a. As discussed in the beginning 
of this subsection, fia ^ 1 is required hj t —>■ fi'j constraints, which implies that 
there are no solutions with 10~^ eV^ < Am^ < 8 x 10~^ eV^, (m^)max < 10 eV^, 
and I tan 2611 > 2.3. 

In summary, we have analysed in detail the issue of obtaining large mixing and 
10~^ eV^ < Am^ < 8 x 10~^ eV^ in the fi — t sector. For all configurations we 
discussed, this is only achieved by carefully choosing the value of specific distances. 
Quantitatively, we find that it is necessary that the value of two a priori independent 
distances be equal to, at least, 1 part in 10^ in order to obtain large mixing in the 
atmospheric neutrino sector. Cases where maximum mixing is "naturally" obtained 
by placing fields "on top" of one another (Cases C and D), were ruled out by flavour 
changing neutral current constraints combined with obtaining the observed values 
for the corresponding charged lepton masses. 

We conclude with a comment on the order one higher dimensional Yukawa 
couplings which we set to unity throughout. Firstly, note that these coefficients 
can be readily "absorbed" by slightly modifying the distance between the lepton 
fields. Explicitly, 



ae-^^' = e-^('')^ (3.31) 
2\na-{ixlf = -(/i/')^ (3.32) 

where ///' is how much the distance has to be changed in order to "absorb" the 
coefficient a. For neutrino masses, (/i/')^ ^ 50, such that, even if a were as small 
as 0.1, \l/l' — 1| ^ 0.05. Secondly, and most importantly, the presence of the order 
one coefficients will not affect our fine-tuning discussions - parameters just have to 

^One may worry what happens if (3 is very large. In this case D2 is tiny, which efficiently 
suppresses this term. 
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be slightly "redefined." For example, in Case A, instead of having to "tune" li to 
I2, one would tune h to I2 plus a small, constant, coefficient dependent on the ratio 
of the higher dimensional Yukawa couplings. The "amount" of fine-tuning remains 
the same. 



3.2 Three Neutrino Families 

In this subsection, we discuss the case of three lepton families, concentrating on ob- 
taining appropriate neutrino mass-squared differences and mixing angles for solving 
the atmospheric and solar neutrino puzzles, while satisfying the constraints from 
the CHOOZ and Palo Verde experiments JTB[.p] 



Unlike the two-family case, in the case of three families it is not simple to perform 
a search for all configurations which potentially yield the appropriate values for 
the oscillation parameters. Furthermore, after some configurations are identified, 
a general analytic discussion of "fine-tuning" among different parameters is not 
illuminating. For these reasons, we decided to approach the situation in a different 
way. 

Initially, we numerically scanned the "location space" of the left-handed lepton 
doublets and the right-handed neutrino singlets, restricting ourselves to branes which 
are fat in only one extra dimension. During the scan, we looked for solutions which 
yield neutrino masses and mixing angles such that all the neutrino experimental 
data could be accounted for by neutrino oscillations [Q. Details of the numerical 
scan are described in the Appendix, and some results are the following: i- As 
expected, obtaining reasonable values for the mixing angles is quite hard, more 
so than the correct order of magnitude for the mass-squared differences, ii- In all 
of the acceptable configurations, the /x-r system falls into one of the configurations 
(see Fig. 0) discussed in detail in the previous subsection, with the same fine-tuning 
characteristics, in- Only configurations which yield hierarchical neutrino masses 
were found, iv- For a significant fraction of the points (~50%) flavour changing 
effects in the charged lepton sector are absent, and the charged lepton mass matrix 
is virtually diagonal. In summary, it was possible to find several configurations which 
yield appropriate neutrino masses and mixing angles for explaining all neutrino data. 

Four different configurations are depicted in Fig. ^.Q These are labelled PLMA, 
PSMA, PLOW and PVAC, alluding to the fact that they will fall into the LMA, 
SMA, LOW and VAC solutions to the solar neutrino puzzle, respectively. Table § 
contains the exact distances for the four configurations.]^ Note that PLMA and 
PVAC have the second and third generation configuration depicted in Fig. ^B), 
while PSMA and PLOW present that depicted in Fig. |](A). 

In order to try to quantify if /how these solutions are finely-tuned, we study the 

"in this paper, we will not try to accommodate the LSND data which is still to be confirmed 
by another experiment, and will only consider three active neutrino species. 

**For all four examples, the brane seems to be too fat. This is a generic feature of almost all 
solutions we find, as discussed in the Appendix, and should be interpreted as a motivation for 
introducing more than one "fat" extra-dimension. 

''"'^see the Appendix for a definition of the "distance parameters." 
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PLMA , , , , , 



PSMA , , , , 

^ 7.7fx ^ 7.7fx J 8.0|a ^ 8.5)J, ^ 

1^. n^n, Ix le 



PLOW , , , , 



le V Ix 



le n^n, Ix \ V 



Figure 4: Four examples which yield neutrino masses and mixing angles which satisfy 
all neutrino data. 



stability of the four configurations depicted in Fig. | by probing "around" them. 
Explicitly, we fix the "smallest" parameter ( I2 in the case of PSMA and PLOW and 
Z3 in the case of PSMA and PLOW) and vary the other ones {li,x,y, and I2 or ^3) 
from Pi — 0.1 to Pi + 0.1 in steps of 0.02, where Pi is the original value of each of the 
four parameters (for a total of 11^ points). The size of the interval is guided by the 
results of our two-generation studies in the previous subsection, i.e., we choose an 
interval which is small enough so that the large atmospheric angle is not disturbed 
too much. 

Fig. § (Fig. I) depicts the resuh of this stability study around PLMA (PSMA). In 
the figures, we plot the values obtained for all the oscillation parameters: Aml^^ x 
tan^6'atm (top left), Am^^^^ x |f/e3p (top right), and Am| x tan^ ^0 (bottom). The 
thin [blue] boxes mark the "allowed regions," as defined in the Appendix, while the 
solid, dark dot indicates the results obtained for PLMA (PSMA). When perturbing 
about PLOW and PVAC, similar results are obtained (indeed, we have repeated 
this procedure for a large number of points). 

The results contained in the figure are quite interesting. First, it should be noted 
that, as expected, the atmospheric parameters do not vary substantially. Numer- 
ically, 56% (54%) of the points remain inside the "atmospheric box", 10~^ eV^ < 
Aml,^ < 10-2 eV2, 1/3 < tan^ O^t^ < 3 in the case of PLMA (PSMA). On the other 
hand, some of the solar parameters and |f/e3p vary drastically. For example, while 
perturbing around PLMA (PSMA), tan^ 6'© varies by 7 (9) orders of magnitude, 
and similar results apply to |f/e3p. Moreover, Am-Q, while somewhat stable around 
PLMA (it varies by 2 orders of magnitude), spans 5 orders of magnitude around 
PSMA. 
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I I I I I I I I I I I I I I 

-3-2-1 2 3 4 

10 10 10 1 10 10 10 10 

Figure 5: Values of Aml^^ x tan^ 6'atm (top, left), Aml^^ x |?7e3p (top right), and 
ArriQ X tan^ ^0 (bottom) which are obtained when perturbing around PLMA (see 
text for details). The results obtained at PLMA is indicated by a solid, dark dot. 
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Figure 6: Same as Fig. |, when perturbing around PSMA. 
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Table 4: Points which yield neutrino masses and mixing angles such that all 
neutrino data can be accommodated. The different points fall (PLMA, PSMA, 
PLOW, PVAC) in different solutions to the solar neutrino puzzle, and are labelled 
accordingly. For the definition of the parameters li,l2,h,x,y, see Fig. || in the 
Appendix. 



configuration 


fill 






fix 




PLMA 


8.0 


-7.7 


7.9 


-16.0 


-15.6 


PSMA 


8.5 


7.7 


7.7 


15.7 


0.0 


PLOW 


-8.4 


7.7 


7.7 


-16.0 


0.0 


PVAC 


-8.8 


-7.7 


8.7 


0.0 


-16.4 



Does this observed numerical instability indicate a fine-tuning among different 
parameters? After all, what we would like to learn from these "stability studies" 
is whether, in order to guarantee that the solar parameters and |?7e3p are inside 
a certain range, the distances between some lepton fields have to be chosen very 
carefully, and whether slight modifications severely spoil the result. In order to 
verify that this is the case, we look at PLMA and PSMA in more detail. The 



lead to the following neutrino Dirac mass 



PLMA-like configurations (see Fig. 
matrix: 

= 



where a, a', b, b', c are related, respectively, to the distances \le—ne\, 
I^T — n^\, \lr — Url- Eq. (|3.34|) leads to 




\le-nr\ 



(3.34) 



/a2 + a'2 



rriym 



a'c 



V 



a'c 





62 bb' 
bb' 6'2 



(3.35) 



From the discussion in the previous subsection, we know that large mixing in the 
atmospheric sector is obtained for 6 ~ 6' > c (Case B). It is also easy to see that, in 
order to obtain Amg < Am^^j^, a, a', c <b, b' . In hght of this, we multiply Eq. ( p.35|) 
by a rotation about the "e-axis" such that sin d = b'/VWTb^ and obtain 



fl 



VO 











COS 9 sin 9 
— sin 9 cos 9 

fa^ + a'^ 
a'c sin 9 








a'c 



(\ 



bb' 







cos^ 
y sin 6* 



62 
bb' 

a'c sin 9 a'c cos 9 

6^ + 6'2 + sin^ 9 c^ sin 9 cos 9 

c^ cos^ 9 





— sin 9 

cos 9 



\ a'c cos 9 (? sin 9 cos 9 



(3.36) 



If \? ^ ^ c^, a'c, it is easy to see that 9 ~ 6'atm (corrections to this lead to a nonzero 
?7e3, for example), as expected. Furthermore, the two other eigenstates are "light" 
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fas lone; "small") and the solar angle is easily calculable: 

c2cos2e-(a2 + a'2) ^ ^ 

cos 2^0 - ^ ^ ' (3.37) 

'{a? + a'2 - c2 cos^ 6Y + 4a'2c2 cos^ 9 



In order to obtain, for example, large mixing in the solar sector, it is necessary 
to satisfy + a'^ ~ cos^ 6* ~ c^/2. This is what happens at PLMA: a = a' = 
g-o. 5x8.0 ^^^^ g _ g-o. 5x7.9 ^ g^ch that a/c 0.45. Therefore, a large solar angle 
is a consequence of finely tuning the distances |/e — ne,r| and 1/,- — n^|. It should 
be noted that the results of the scan involve further complications, such as c ~ 6', 
which leads to large values of Am^ and significant corrections to the atmospheric 
angle. 

The PSMA-like configurations (see Fig. ^) lead to the following neutrino Dirac 
mass matrix: 

/e 

rriy = 



b b' \, (3.38) 
a b" b'" 



where e,a,b,b',b",b"' are related, respectively, to the distances |/e — ne|, \It — ""-el, 
l^n — "^mI' I^m ~ — \lr — rirl- From Fig. ^, it is easy to see that e <^ a <^ 

b, b', b", b'". Eq. diH) leads to 

/ ae 

m^ml ~ b'^ + b'^ bb" + b'b'" \ . (3.39) 
\ae bb" + b'b"' b"^ + b"'^ + a^ 

As before, the atmospheric angle and mass difference result from diagonalizing the 
fj,-T sector. As we did in the PLMA case above, we can "undo" the atmospheric 
rotation and obtain 



I e^ aesmO aecosO 

aesinO ml + a"^ sin^ 9 a'^ sin 9 cos 6 I, (3.40) 
Uecos^ Jsin^cos^ m^ + a^cos^^, 



where m\,m\ are the lightest and heaviest eigenstates obtained after diagonalizing 
the ii-T sector. They are given by 



'2 ^ 2 



Large mixing in the atmospheric sector requires b ^ b' ^ b" b'" (see Case A in the 
previous subsection), and, in the case of hierarchical neutrinos (the case of interest 
here), Am^^-j^^ ~ > mf. In this case, it is trivial to compute the solar angle and 
mass squared difference, keeping in mind that a^ ^ e^ and that cos^ Oatm — 1/2: 

Am| ~ a^/2 + m\, (3.42) 



tan^0~sin^Q ~ ^'^l^ ^ . (3.43) 
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Exactly at PSMA, ml = 0, and AttIq and ta.n^ 6q are governed by a and e. As 
long as this is the case {ml <C a^), varying li, I3 and x (which changes the distances 
\le — Uel and \It — Uel) only perturbs around PLMA. However, when mf grows (this 
happens for y ^ 0) the situation changes dramatically, and Am^ grows, while 
tan^ 6'q decreases sharply. This explains the "cloud" of points in Fig. |^(bottom) at 
10-^ eV^ < Ami ~ 3 x lO'^ eV^, 10^^° < tan^^o < 10^^^ In summary, the solar 
parameters depend very strongly on y, the distance between and rir- More so 
than the atmospheric parameters! Note that a very similar behaviour is observed 
when perturbing around PLOW, where the /x-r fields also fall into the configuration 
discussed in Case A in the previous subsection. 

We conclude this subsection with a summary of the results and a word of caution. 
We have numerically scanned the one- dimensional parameter space of three neutrino 
generations and encountered, out of many "good" candidates (which yield small 
neutrino masses), only a handful of points which could potentially accommodate the 
neutrino data. As expected from our qualitative discussion at the beginning of this 
section, the great majority of the points (> 99.9%) fails to yield large enough mixing 
in the solar and/or atmospheric sectors. The situation is particularly constrained 
after |f/e3p is required to be small. All the points identified had configurations and 
fine-tuning characteristics similar to the ones discussed in the previous subsection 
(see Fig. |^), with different choices for the position of the first generation fields. 

We selected four examples consistent with the different solutions to the solar 
neutrino puzzle, and studied the stability of these solutions in order to address 
how "special" the distances between lepton fields have to be in order to obtain a 
particular result. We find that small modifications (of order a few percent of the 
typical distance scale) seem to maintain the large mixing in the atmospheric sector, 
in agreement with the discussion in the previous subsection, while changing the solar 
parameters and |?7e3p by many orders of magnitude. For two specific cases we iden- 
tified the origin of this effect. Therefore, it seems that obtaining appropriate solar 
parameters requires tuning more parameters (or the same parameters more tightly) 
than what is already required in order to obtain large mixing in the atmospheric 
sector. 

It should be emphasised that we have not attempted a detailed study of the three 
by three neutrino mass matrix in order to study correlations between the various 
distances, but decided to base our conclusions on a few "case studies." It is not 
clear to us whether such a study is practical and/or illuminating. Therefore, the 
conclusions we reached, although quite plausible and reasonable, should be read 
with some caution. 

4 Neutrino Masses without Right-Handed Neu- 
trinos 

As mentioned in Sec. 2, the traditional see-saw mechanism is not an option for 
generating very small neutrino masses in the case of models with a small quantum 
gravity scale. Nonetheless, one may still obtain small Majorana neutrino masses in 
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the ADD scenario. The idea, proposed in |T8[ , is to take advantage of the "infrared 



desert" [|r^ to obtain very small symmetry breaking effects. More explicitly, the idea 
is to impose lepton number as a conserved global symmetry, which is broken at some 
far away brane. Lepton number violating effective operators would be generated in 
our brane, suppressed by e""*^ where m is some typical mass scale and r is the 



distance between the branes (see [|T8| , p!9[] for details). 



For our purposes, it is enough to assume that the following Majorana neutrino 
mass matrix is generated after electroweak symmetry breaking: 



(4.1) 



where mMaj is the overall Majorana mass scale, ap^ are the higher dimensional order 
one couplings.]^ The distances |/j — lj\ are the same distances between the different 
left-handed lepton fields discussed in the previous sections. 

The issue we would like to deal with in this section is whether phenomenologically 
acceptable Majorana mass matrices can be obtained, if one assumes that the absence 
of flavour violating decays of the muon and the tau and the hierarchy of the charged 
lepton masses are explained by the AS scenario. 

In the absence of the flavour changing constraints, it is clear that solutions can 
be found: if one places all lepton doublets on top of one another, large mixing 
in the solar and atmospheric sector arises quite naturally, and even a small \Uez\^ 
and moderate hierarchy in the masses-squared can be obtained With the 



flavour changing constraints, however, the situation is more involved, since the lepton 
doublet fields cannot sit on top of one another, as we argued in Sec. 3. Note that 
if all distances were |/a — //j] ~ 2/i~^ all off-diagonal terms would be one order 
of magnitude smaller than the diagonal terms,Q and obtaining large mixing in the 
atmospheric sector, for example, would require some fine-tuning among the diagonal 
elements, which we do not consider acceptable. Note that the situation deteriorates 
very quickly as — Ipl grows, due to the Gaussian dependency of the off-diagonal 
coefficients. Therefore, the first issue we will address is how small can we make 
|/o — without running into contradictions with flavour violating charged lepton 
decays. 

From the previous sections, we learn that \le — ~ 3/i~^, while \lr — le,^l\ ~ 
Note that the lower limits cannot be simultaneously saturated due to the 
(much more stringent) /i — 67 bounds and the requirement that the charged lepton 
mass hierarchy is explained within the AS scenario. For example, if we require 
\It — ~ in one extra dimension, — /e| > 7/i~^, while \le — ~ 9>ii~^. 

*We are assuming, according to the general philosophy of the rest of the paper, that no flavour 
structure comes from the lepton number violating sector. Such a structure (conserved flavour 
symmetries, for example) would manifest itself by imposing relations between the ap^. 

^Indeed, in this case large mixing would also come from the charged lepton sector. 
According to the authors of jl^, gauge interactions may significantly alter this picture. For 
simplicity, however, we will ignore these effects, which seem to be very dependent on the fermion 
localising mechanism. 
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In two or more extra dimensions, it is possible to obtain 1/,- — le,\ ~ 2.5/i~^ while 
keeping \le — lfj,\ ~ 3/i~^ (in which case, \le — Ifil is very close to its lower bound). It 
should be kept in mind that this situation is very "special," and the branching ratios 
for fj, — > eee and r {e,fj,)j are very close to the current experimental bounds p. 

Therefore, the "best" possible scenario seems to be: the e/x-element is suppressed 
by at least e~°'^^^ ~ 0.01 with respect to the diagonal elements, while one of the 
other off-diagonal elements (we choose it to be the /xr-element) can be suppressed by 
only e~°-^^^ ~ 0.6 with respect to the diagonal terms. The remaining off-diagonal 
element (we will choose it to be the er-element) is suppressed by 0.05. At this point, 
one may ask whether it is possible to find, without carefully choosing the order one 
af3^ coefficients, phenomenologically viable neutrino masses and mixing angles. 

In order to answer this question, we perform the following "search." We parametrise 
the neutrino Majorana mass matrix by 



= mMaj 



a 



ee 



V 




a^f, a^r X 0.5 , (4.2) 



where we randomly choose values for Ci(^ee),{^in),iTT),{pT) between 0.3 and 3.0,|^ while 
randomly choosing values for e between and 0.1. The e/i-entries have been set 
to zero for simplicity. We have verified that the results we obtain are not changed 
in any significant way if e/x-elements of order 0.01 are also included. Note that 
we still have one free parameter, namely mMaj- The quantity mMaj ■ «ee is directly 
constrained by searches for neutrinoless double beta decay, which currently requires 
\Ml'^\ < 0.35 eV 0, |20[. In our analysis, we will circumvent the dependency on mMaj 



by computing the ratio of the mass-squared differences 

Aw? 

''atm 



ratio = , ":® , (4.3) 



where Am^'s and mixing matrix elements are defined as in Eqs. ( [A.l| ) (see also the 
remarks which follow the equations for a comment on inverted hierarchies). 

Figs. 0(top) depict the results obtained for 10001 randomly selected points. 
Fig. ^(top,left) depicts the obtained values of ratioxtar? 6'atm, while Fig. ^(top, right) 
depicts the obtained values of |C/e3p x tan^6'0. As one could have expected, the 



atmospheric angle is generically large [p!3| , |TJ] even though the /xr-element of Eq. ( ^^21 



is typically suppressed by a factor of two with respect to the diagonal elements, and 
ratio is indeed usually of order one. Furthermore, |t/e3p and tan^ ^0 are usually 
small, even though large values of tan^ 6q can be obtained.^] 

More significant than the figures, perhaps, are the following numerics: 51% of the 
points yield 1/3 < tan^ ^^atm < 3, while 22% (1.5%) of the points yield ratio < 10^^ 

§we are neglecting, for simplicity, possibly large phases associated with the order one higher 
dimensional Yukawa couplings Q . 

^Note that there are some points for which |J7e3p — 1- These are expected, and correspond 
to cases when aee is bigger than the other coefRcients. In this case, the largest eigenvalue is 
proportional to aee and has a very large electron- type component. All these points are incompatible 
with the current neutrino data, of course. 
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Figure 7: Neutrino mixing parameters obtained for 10001 randomly generated 
Majorana mass matrices, imposing constraints from the absence of flavour violating 
charged lepton decays (see text) displayed in the following manner: ratio = 
AttIq/ Aml^^ X tan^ 6'atm(^e/t) and |?7e3p x tan^ 9q (right). In the top, the branching 
ratios for flavour changing lepton decays are close to the current experimental bounds 
0, while the bottom figures explore the behaviour of the oscillation parameters as 
the prediction for the branching ratios for flavour changing charged lepton decays 
decreases (see text for details). 
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(10~^), and 38% of the points yield 10~^ < tan| < 1 (note that the values of ratio 
we obtain are not small enough to accommodate the LOW of VAC solutions to the 
solar neutrino puzzle). Combining all constraints, 1.5% satisfy 1/3 < tan^ 6'atm < 3, 
10"^ < ratio < 10-^ 0.1 < tan| < 1, |[/e3p < 0.1 (LMA-like solution) and 0.1% 
satisfy 1/3 < tan^ ^^tm < 3, lO'^ < ratio < 10'\ IQ-^ < tan| < 2 x 10-^ 
|C/e3p < 0.1 (SMA-like solution). 

What do these numbers imply? While most points do not yield phenomenologi- 
cally acceptable parameters, there is a slightly more than 1% chance for obtaining 
the LMA solution with completely random order one Yukawa couplings. We do 
not consider this too "unlikely," but leave definitive conclusions to the reader. 
One concrete statement, however, is that obtaining the SMA solution under these 
conditions is certainly less likely (15 times). The reason for this is that the SMA 
solution requires a larger hierarchy for the mass-squared differences, which is very 
hard to obtain with order one random Yukawa couplings. 

What happens as the separation between left-handed lepton fields increases? 
Fig. ^(bottom, left) depicts the obtained values of ratio x tan^ 6'atm if the same 
procedure is repeated, except that the suppression factor for the /ir-element is 
increased by a factor of 10 (from 0.5 to 0.05). The results are rather striking. 
Unlike the previous case, depicted in Fig |^(top, left), large atmospheric mixing 
becomes more "rare" (11% of the points yield 1/3 < tan^ 6'atm < 3), especially 
when ratio is required to be small (0.04% percent (none) of the cases satisfy 1/3 < 
tan^ 6'atm < 3 and ratio < 0.1 (0.01)!). Therefore, in light of the observed hierarchy 
of neutrino mass-squared differences, large mixing in the atmospheric sector becomes 
exceedingly unlikely as the separation between and increases. 

On the other hand. Fig. |^(bottom, right) depicts the obtained values of |t/e3p x 
tan^ 6'0 when the original procedure is again performed, but restraining e < 0.01. 
When compared to Fig. [^(top, right), it is easy to see that large mixing in the 
solar sector happens for a smaller fraction of all randomly generated matrices. The 
number of potentially good candidates is still reasonable: 28% of all matrices yield 
10^^ < taug < 1. Here, however, the probability of obtaining LMA-like solutions 
(see discussion related to Fig. 0(top)) decreases (0.3% compared to 1.5%), while 
the number of SMA-like solutions remains the same. Of course, this behaviour 
is expected, and the situation deteriorates very rapidly as the upper bound on e 
decreases. 

In summary, if neutrinos have a naturally small Majorana mass matrix, the AS 
solution to flavour changing charged lepton decays and the hierarchy of charged lep- 
ton masses naturally predicts that the neutrino masses are not strongly hierarchical 
and that mixing angles are rather small. In order to obtain large enough atmospheric 
mixing angles, we were forced to live very close to the current experimental bounds 
imposed by the searches for charged lepton flavour violation. We dare not make 
any precise predictions given the qualitative nature of our estimates, but if flavour 
violating fi and (especially) r decays are not observed by the next generation of 
experiments [^, the scenario discussed in this section will be severely constrained, 
if not entirely ruled out. 

This is to be contrasted to the study in the previous section, where neutrino 
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mixing angles were also gcncrically predicted to be small, while the masses were 
hierarchical. In order to obtain large mixing angles, we were forced to impose 
specific positions for the different leptonic fields. All solutions we found were either 
immediately ruled out by searches for lepton flavour violating processes, or predicted 
that their branching ratios were severely suppressed, out of reach of any foreseeable 
experiment. The price paid for these solutions was having to finely tune a priori 
unrelated distances to more than 1 part in 10^. 

5 Summary and Conclusions 

We have analysed the issue of charged lepton fiavour violation and neutrino masses 
in the fat-brane paradigm proposed by Arkani-Hamed and Schmaltz (AS scenario). 
The AS scenario is particularly suited for explaining the absence of lepton flavour 
violating muon and tau decays, and can also naturally explain why neutrinos are 
more than ten orders of magnitude lighter than the top quark, if right-handed 
neutrinos are introduced in the brane. 

We found that, indeed, the branching ratios for flavour violating muon and 
tau decays can be very easily suppressed to levels below the current experimental 
bounds, and that very small neutrino masses can be obtained for separations of 
order 10/i~^. The AS scenario, however, seems to have a harder time accounting 
for the large neutrino mixing which has been observed in the atmospheric neutrino 
data. This is a peculiar feature of the AS scenario for fermion masses: it naturally 
accommodates large mass hierarchies and small mixing angles, while it seems to 
require additional "structure" {e.g., different pairs of fields have to be separated by 
very similar distances) in order to explain large mixing angles. 

We studied the case of two lepton families in detail. The situation seems "finely 
tuned" when the the bounds imposed by charged lepton flavour violation are also 
taken into account. We showed that the charged lepton mass matrix is constrained 
to be quasi- diagonal, and that, in order to obtain large atmospheric mixing and 
the appropriate Am^, distances which are a priori unrelated are required to agree 
with one another better than one part in 100. We proceeded to analyse the case 
of three lepton flavours, and while we found that all the different solutions to the 
solar neutrino puzzle can be accommodated while keeping the absolute value of the 
Ue3 element of the leptonic mixing matrix small, wc argued that the amount of fine- 
tuning required to simultaneously satisfy all neutrino data is "larger" than in the 
two-family case. 

We then proceeded to discuss the effect of explaining the absence of flavour 
changing tau and muon decays and charged lepton masses in the AS scenario on 
theories where small Majorana neutrino masses arc generated by breaking lepton 
number in a far away brane. In this case, in order to obtain large mixing in the 
atmospheric sector, we were forced to "approximate" different lepton doublet fields 
to the point that the current experimental upper bounds on the branching ratios of 
some rare tau and muon decays were almost saturated. Therefore, we expect these 
rare decays to be observed in the next round of experiments, if such a scenario were 
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indeed realised in nature. Furthermore, very hierarchical neutrino mass-squared 
differences were not attainable, meaning that the solar neutrino puzzle would have 
to be solved either by the SMA or the LMA solutions. 

In conclusion, we would like to emphasise that, in spite of the fine-tuning prob- 
lems we highlighted in Sec. 3, the AS scenario should be regarded as a novel 
mechanism for understanding small mixing angles, which does not necessarily require 
the existence of fat branes or even large extra dimensions. Similar results are 
obtained in models with many different "thin" branes which are separated in the 
extra dimensions. Furthermore, if the extra dimensions are not large (and the 
hierarchy problem is solved by some other means), our results for neutrino masses 
and mixing angles would still apply, except for the flavour changing constraints 
from charged lepton decays, which would be absent. Finally, the fine-tuning which 
we pointed out could (and should, in our opinion) be interpreted as a challenge to 
be faced by candidates for the localising mechanism and/or the structure of the 
compact dimensions. 
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A Numerical Scan of the 3 Generation Parameter 
Space 

In this appendix, we describe a numerical scan over the three generation, one extra 
dimension location space of the lepton doublets and right-handed neutrino singlets. 
In this case, specifying five distances suffices to determine the relative distance 
between the six different fermion fields. These distances are pictorially defined in 
Fig. and are: /i, I2, the distances between la and the distance 

between and rig; y, the distance between rij. and n^. Note that the sign of four 
out of the five distances is also meaningful (for example, when y < 0, is "to the 
right" of rir). 

In the scan, we varied h {I2) from — 20/i^^ to 20^^^ (— 20.1/i^^ to 20.1/i^^ ) in 
steps of 0.2, X and y from —20^^^ to 20^~^ in steps of 0.4, and from 0.1^^^ to 
20.1/i~^ in steps of 0.2. Out of all the points, we only consider these which yield 
neutrino mass matrices such that {m^ml)ij < 10 eV^ for all i,j = e,/i,r. This is 
done in order to guarantee that all neutrino masses squared are roughly less than 
10 eV^. One consequence of this constraint is that for all U, ^ 7.0, while certain 
bands for x and y are excluded. 
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Figure 8: General three family configuration in one extra dimension. There are five 
independent distances, which are chosen to be li,l2,h,x,y. Note that in order to 
span the entire location space, four out of the five parameters should take positive 
and negative values, while the sign of the fifth one remains fixed. 

Of the remaining matrices, we identify the ones whose eigenvalues and eigenvec- 
tors fall in the following window: 

\UeZ? < 0.1, 

l/3<(|?7^3|/|f/r3|)'=tan2^^,tm<3, 

lO-5<(|t/e2|/|f/ei|)' = tan2^0<lO, 
10-3 eV^ < {ml - ml) = Am^^^ < lO'^ eV^, 
10-^° eV2 < {ml - ml) = Ami < 10"^ eV^, (A.l) 

where we define the neutrino mixing matrix as Ui = Uaii^a, where z/q,, a = e,fi,T 
are weak eigenstates and z/j, i = 1,2,3 are mass eigenstates, with masses mi. The 
eigenstates are organised in ascending order of mass-squared. When m| — ml > 
ml — ml (normal hierarchy), the window defined above contains the region of the 
parameter space which satisfies all of the experimental constraints . If m^ — ml < 
ml — ml (inverted hierarchy) the same window can be used after relabelling the 
mass eigenstates 1 ^ 3 — 2 — » 1, and defining {ml — m|) = Am^^^. It should be 
noted that during the numerical scan we did not find any matrix which yields an 
inverted mass hierarchy and is consistent with the constraints imposed above. 

The results of the scan are best described in words. Out of roughly 2 x 10^ initially 
selected matrices, only 440 fall within the window defined above. The events are 
rejected at the following rate: 99% fail the mass constraints. From the surviving 
points, 96% fail to satisfy the tan^ 6q constraint. From the remaining points, 97% 
fail the tan^ 6'atm constraint, while 98% of the remainder fall out of the |t/e3p bound. 
Note that only one out of roughly 5 million events survives all the constraints we 
define above. 

Many comments are in order. First, in should be readily noted that the two 
"angle constraints" combine to remove 99.9% of the candidate matrices, ten times 
more than the "mass constraints." Second, it should be pointed out that the mass 
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constraints and the angle constraints "commute", in the sense that, independent of 
the order in which they are apphed, the percentage of the points which is removed 
is the same. This imphes that there is no correlation between the values of the 
masses and the values of the mixing angles. The same is not true for the tan^ 
and tan^ 6'atm constraints. Third, the |f/e3p constraint is highly correlated with the 
angle constraints. If applied first, the |f/e3p constraint only cuts a small fraction 
of all the matrices, while if applied last, it removes a very significant fraction of 
the events which pass the angle constraints. Therefore, events which have a small 
{Uesl generically have very small mixing in the solar and/or the atmospheric sector. 
This can be understood from the qualitative discussion included in the beginning 
of Sec. 3: the AS scenario prefers small mixing angles, and a generic set of points 
will yield a mixing matrix which has only "Os" and "Is" (one 1 for each column 
and row). Therefore, if we impose |t/e3p small, the generic matrices which are left 
have some other Uas and Uei which is close to one, which necessarily means that 
tan^ ^0,atm are either very large or very small. 

The selected points yield a varied spectrum of solar angles and masses-squared, 
and no value of Am^ and tan^ ^0 is strongly preferred. All points, without exception, 
choose the second and third generation fields to fall into either the configuration 
depicted in Fig. ^(A) or Fig. 0(B) (note that the configurations where and Ur 
are exchanged are, of course, also present in equal numbers), while li and ni are 
placed in various locations. Furthermore, the choices for I and 6 in the case of the 
configuration Fig. |^(A), and /21 and I22 in the case of configuration Fig. §(B) are 
exactly the ones obtained in the last line of Tables |l| and 0, respectively. This implies, 
as one would naively expect, that in order to obtain large mixing in the "atmospheric 
sector" the same tuning of parameters as the ones discussed in Sec. 3.1 is required. 

Fig. P(left) is a scatter plot of the values of I2 and li of the points which pass 
all of the constraints mentioned above. The light [green] points are the ones which 
also pass the constraints that two lepton doublets are not too close to one another 
(flavour changing constraints): 1^,- — le,fj.\ ^ 1/^""^ and — le\ ^ (see Fig. |^). 
This constraint is sufficient for satisfying all flavour changing constraints described 
in Sec. 2 and also guarantee that the mixing in the charged lepton sector is negligible. 
218 out of the 440 points satisfy this extra requirement. 

Fig. ^( right) is a histogram of the largest distance between two lepton flelds, and 
can be interpreted as a lower bound on the brane "thickness." The solid painted 
histogram is for points which satisfy the flavour changing constraints. Note that 
only "very thick" branes are obtained, especially for matrices which do not have 
two Is on top of each other. This could be the subject of serious concern, since 
branes should not be too "fat" if one wants to solve the gauge hierarchy problem 
and still have the gauge bosons propagate freely in the entire fat-brane (see for 
details). We prefer, instead, to view this as an indication that one extra dimension 
seems to be problematic, and that if one considers branes which are fat in more 
dimensions similar conflgurations can be obtained for much "thiner" branes. This 
is very easy to see. All one has to do is rotate a number of lepton flelds, keeping the 
"smaller" distances flxed, while making sure that the "larger" distances remain large. 
Because of the Gaussians dependency of the Yukawa couplings on the distances, 
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Figure 9: Scatter plot of I2 x li of all the 440 points which satisfy all constraints 
imposed by the neutrino data (left, see text). The light [green] points also satisfy 
the "flavour changing constraints," |/r — ^e,/i| > 1/^"^ and \lf^ — le\ > 3fi~^. Histogram 
of the brane-thickness for the same 440 points (right). The painted histogram 
corresponds to events which satisfy the flavour changing constraints. 

these configurations will be, in practice, indistinguishable from the one-dimensional 
ones we discuss here. 

One may worry whether configurations Fig. ^C) and Fig. |^(D) can yield proper 
solar parameters, since these configurations never show up in our general scan. 
In order to check this, we numerically searched for solutions by starting from 
configurations Fig. |^(C) and Fig. |^(D), choosing distances such that the atmospheric 
neutrino puzzled is solved and varying the parameters Zi and x (see Fig. Indeed, 
we are able to find a large number of solutions as long as is not on top of Ir, 
but extremely close. This is easy to understand: when is exactly on top of Ir, 
the neutrino mass matrix has a zero eigenvalue, whose eigenvector is very close to 
l/^/2{\u^) — l^r))- This implies that the lUf.) state is contained in the two heaviest 
mass eigenstates, which makes it very hard to solve the solar neutrino puzzle. The 
only alternative is to choose the other distances such that an inverted mass hierarchy 
is obtained, but such a situation is extremely finely-tuned.Q It should be noted 
that all solutions found starting from the configurations H(C) and (D) are ruled 
out by flavour changing constraints. Indeed, as was pointed out in Sec. 3.1, the 
configuration Fig. ^(D) is already rule out in the case of two generations, while 

*One may wonder whether an SMA-hkc solution to the solar neutrino puzzle cannot be obtained 
by varying the parameters slightly. This is not the case since i^e is going to be predominately heavy 
and even though sin^ 26'© ^ 1 this point belongs to the "dark side" of the parameter space, 
tan^ 6*0 > 1. 
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satisfactory solutions with the configuration Fig. ^(C) are extremely finely tuned. 

Another worry is whether much thinner branes cannot be obtained. One pos- 
sibility, for example, is to have all right-handed neutrinos very close to each other. 
Since we are restraining ourselves to only one dimension, these configurations are 
very hard to obtain, and the fiavour changing constraints make them even harder to 
find. Nonetheless, we have been able to find (numerically) examples which fit this 
description and do not yield too large charged lepton mixing or fiavour changing 
effects. 
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